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The standard mean-field theory for the Mott insulator-superfluid phase transition is not sufficient to describe 
the Mott insulator-paired superfluid phase transition. Therefore, by restricting the two-species Bose-Hubbard 
Hamiltonian to the subspace of paired particles, and using perturbation theory, here we derive an analytic mean- 
field expression for the Mott insulator-paired superfluid transition boundary. 
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I. INTRODUCTION 

Following the recent observation of Mott insulator- 
superfluid transition with ultracold atomic Bose gases loaded 
into optical lattices tlj-lj], there has been an intense theoretical 
activity in analyzing many Hubbard-type lattice models flU]. 
Among them the two-species Bose-Hubbard model, which 
can be studied with two-component Bose gases loaded into 
optical lattices, is one of the most popular. This is because, 
in addition to the Mott insulator and single-species superfluid 
phases, it has been predicted that this model has at least two 
additional phases: an incompressible super-counter flow and 
a compressible paired superfluid phase lla- lllh . 

Our main interest here is in the latter phase, where a di- 
rect transition from the Mott insulator to the paired super- 
fluid phase (superfluidity of composite bosons, i.e. Bose-Bose 
pairs) has been predicted, when both species have integer fill- 
ings and the interspecies interaction is sufficiently large and 
attractive. In this paper, we derive an analytic mean-field 
expression for the Mott insulator-paired superfluid transition 
boundary in the two-species Bose-Hubbard model. The re- 
maining paper is organized as follows. After introducing the 
model Hamiltonian in Sec. [TTJ first we derived the mean-field 
theory in Sec. Ill Al and then presented typical phase diagrams 
in Sec. MB I A brief summary of our conclusions is given in 
SecHm 



II. TWO-SPECIES BOSE-HUBBARD MODEL 

The two-species Bose-Hubbard Hamiltonian is given by, 

H = - %<A>&7> + -j-niA^i," - l ) 

where the pseudo-spin a = {f, 1} labels the trapped hyperfine 
states of a given species of bosons, or labels different types 
of bosons in a two-species mixture, iy )CT is the tunneling (or 
hopping) matrix between sites i and j, b\ a (bi. a ) is the bo- 
son creation (annihilation) and rii^ = b\ a bi. a is the boson 
number operator at site i, U aa > is the strength of the onsite 



boson-boson interaction between a and a' components, and 
fji(j is the chemical potential. In this manuscript, we consider 
a <i-dimensional hypercubic lattice, for which we assume a 
is a real symmetric matrix with elements ty i<7 = t a > for 
i and j nearest neighbors and otherwise. We take the in- 
traspecies interactions to be repulsive ({U^, E/4,4.} > 0) and 
the interspecies interaction to be attractive (Ufi < 0) such 
that U^Un > Ul,, to guarantee the stability of the mixture 
against collapse. 

For sufficiently attractive f/fj., it is well-established that |0- 
[TTll instead of a direct transition from the Mott insulator to 
a single particle superfluid phase, the transition is from the 
Mott insulator to a paired superfluid phase (superfluidity of 
composite bosons, i.e. Bose-Bose pairs). In fact, in the limit 
when {t^, t^} — > 0, it can be shown that the transition is from 
the Mott insulator to a paired superfluid phase for all U^i < 
0|12|. 



A. Mean-field theory 

In the single-species Bose-Hubbard model, the standard 
mean-field theory, where the boson creation and annihilation 
operators are approximated by their expectation values, e.g. 
bi.a = (bi,a) + 5bi, ai has proved to be very useful in under- 
standing the qualitative features of the Mott insulator-single 
species superfluid phase transition |5|]. This is simply be- 
cause the transition is driven by the first-order hopping ef- 
fects. However, the Mott insulator-paired superfluid transi- 
tion is driven by the second-order hopping effects, and there- 
fore, the standard mean-field theory is not sufficient. This 
difficulty could be overcome by restricting the Hamiltonian 
to the subspace of paired particles, and including the second- 
order hopping effects through second-order perturbation the- 
ory iisuJ. 

For the Hamiltonian given in Eq. (Q~|), we have recently 
calculated the two-particle and two-hole excitation energies 
(i.e. energy costs for adding and removing two-particles, 
respectively) up to third order in the hoppings. Assuming 
{U aa , \Ufi\, 2U aa + Ufi} ^> t a , the two-particle excitation 
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energy was found to be 1 1 

Ep = U n (rif + n± + 1) + ^ (Uaafla - Mff) 
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where z — 2d is the coordination number. Similarly, the two- 
hole excitation energy was found to be lfl2tl 

Eh. = -U n (n t + - 1) - [U arT (n a - 1) - 
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The accuracy of Eqs. (O and OJ are checked via exact small- 
cluster (two-site) calculations. In addition, in the limit when 
it = h = *> Un = Un = Uo, Ufi = U', n-j- = nj. = no, 
A*t = A*4- = A 4 ' anc ^ = 2 (or d = 1), Eq. (O is in complete 
agreement with Eq. (3) of Ref. ifloll . providing an independent 
check of the algebra. We note that, unlike the usual Bose- 
Hubbard model where t„ scales as 1/d when d — > oo, here to- 
must scale as 1/ \fd when d — > oo. 

Given the two-particle and two-hole excitation energies, the 
mean-field phase boundary for the Mott insulator-paired su- 
perfluid transition is determined by (see Ref. flUl for a similar 
calculation) 
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where c p = -2(n t + l)(n+ + l)zt^;//7 n and c h = 
— 2n-|-n^zi-f-t^/(7-t-^. We note that, in the limit when tf = tj. = 
J, t/ff = C/44 = C/, J/fi = W ~ —(7, n-j- = ni = m, 
and fi-f — /ij. = /i, Eq. (0| reduces to Eq. (5) of Ref. fPUl 
(after setting Unn = there). However, the terms that are 
proportional to tftf are not included in their definitions of the 
two-particle and two-hole excitation energies. Solving Eq. (0]i 
for fif + hi, we obtain 



Mt + H = 5 



a p - a h ± ^ (a p + a h ) 2 



4c p c h 



(5) 



where a p = £ p + fif + a h = E h - fif - ^4., and ± 
signs correspond to the two-particle and two-hole branches, 
respectively. Equation (O is the mean-field expression for the 
Mott insulator-paired superfluid transition boundary, and it is 
the main result of this paper. 



B. Typical phase diagrams 

In this section, we present typical phase diagrams in the 
fif + hi versus y/ztf plane, obtained directly from Eq. ©. 



Similar to the usual Bose-Hubbard model, as hopping in- 
creases from zero, the range of the chemical potential about 
which the ground state is a Mott insulator decreases, and the 
Mott insulator phase disappears at a critical value of hopping, 
beyond which the system becomes a paired superfluid. 
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FIG. 1: (Color online) The Mott insulator-paired superfluid phase 
transition boundaries are shown for the first lobe, i.e. — 714. = 1, 
when tf = ti- In (a) Ufi = — 0.95f7-f-f, and in (b) Un — Uff. 
Note that the red curves correspond to the same data in both of these 
figures. 

For instance, in Fig.Q] we show the Mott insulator-paired 
superfluid phase transition boundaries for the first lobe, i.e. 
nf = ni = 1, when tf = £4. In Fig. 0I a X where we set 
Ufi = — 0.95Uff, it is clearly seen that decreasing Un favors 
the paired superfluid phase, as intuitively expected. While in 
Fig. |TJb), where we set Un = Uff, it is clearly seen that 
decreasing the strength of Ufi favors the Mott insulator phase 
(see the explanation below). We also note a weak re-entrant 
quantum phase transition in both of the figures. 

In addition, in Fig. [2] we show the Mott insulator-paired 
superfluid phase transition boundaries for the first lobe, i.e. 
rif = ni = 1, when Un = Uff. In Fig. |2ja), where we set 
Ufi = —0.9bUff, it is clearly seen that increasing tj, favors 
the paired superfluid phase, as intuitively expected. While in 
Fig. |2jb), where we set t± = O.Olij,, it is clearly seen that 
decreasing the strength of Ufi again favors the Mott insulator 
phase. However, compared to Figs. [TJa) and QlbX we note 
that the re-entrant quantum phase transition is much stronger 
in these figures. 
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which may indicate an instability toward a collapse. In ad- 
dition, from Eq. [3] we expect a re-entrant quantum phase 
transition when ~(2n t n i /U n )zt t t i - Y,A n l/ u n ~ ( n l ~ 
l)/{2U aa + U n ) + 2n a {n a + 1)/U aa ]ztl < 0, which oc- 
curs beyond a critical Ufp When this expression is nega- 
tive, its value is most negative for the first Mott lobe (i.e. 

= ni — 1), and therefore the effect is strongest there. 
However, its value increases and eventually becomes positive 
as a function of filling, and thus the re-entrant behavior be- 
comes weaker as filling increases, and it eventually disappears 
beyond a critical filling. 
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FIG. 2: (Color online) The Mott insulator-paired superfluid phase 
transition boundaries are shown for the first lobe, i.e. n-j- = nj. = 1, 
when U u = U n . In (a) f/ n = -0.95(7 n , and in (b) t± = 0.01t+. 
Note that the red curves correspond to the same data in both of these 
figures. 



III. CONCLUSIONS 

In this paper, by restricting the two-species Bose-Hubbard 
Hamiltonian to the subspace of paired particles, and using per- 
turbation theory, we derived an analytic mean-field expression 
for the Mott insulator-paired superfluid transition boundary. 
We found that, for small Ufi, the location of the Mott in- 
sulator tip increases as a function of Ufi, because the pres- 
ence of a nonzero is what allowed this state to form in 
the first place. However, when the strength of U^i is larger 
than some critical value (approximately the lo- 

cation of the tip decreases, and it eventually vanishes exactly 
when Ul< — U^U^, which may indicate an instability to- 
ward a collapse. Given that the interspecies interaction can be 
fine tuned in ongoing experiments, e.g. 41 K- 87 Rb [13, 14] or 
homonuclear IT15ll mixtures, via using Feshbach resonances, 
we hope that our predictions could be tested with ultracold 
atomic systems. 



Our results are consistent with the expectation that, for 
small U^i, the location of the Mott insulator tip increases as 
a function of Uf±, because the presence of a nonzero U^i is 
what allowed this state to form in the first place. However, 
when the strength of Ufi is larger than some critical value (ap- 
proximately yjUff Un/2), the location of the tip decreases, 
and it eventually vanishes exactly when U2> = U^Un, 



IV. ACKNOWLEDGMENTS 

The author thanks Christian Trefzger for correspondence. 
This work is financially supported by the Marie Curie In- 
ternational Reintegration Grants (FP7-PEOPLE-IRG-2010- 
268239) and Scientific and Technological Research Council 
of Turkey (TUBITAK). 



[1] M. Greiner, O. Mandel, T. Esslinger, T. W. Hansch, and I. 

Bloch, Nature (London) 415, 39 (2002). 
[2] T. Stoferle, H. Moritz, C. Schori, M. Kohl, and T. Esslinger, 

Phys. Rev. Lett. 92, 130403 (2004). 
[3] I. B. Spielman, W. D. Phillips, and I. V. Porto, Phys. Rev. Lett. 

98, 080404 (2007). 
[4] I. B. Spielman, W. D. Phillips, and I. V. Porto, Phys. Rev. Lett. 

100, 120402 (2008). 
[5] I. Bloch, I. Dalibard, and W. Zwerger, Rev. Mod. Phys. 80, 885 

(2008). 

[6] A. B. Kuklov and B. V. Svistunov, Phys. Rev. Lett. 90, 100401 
(2003). 

[7] E. Altman, W. Hofstetter, E. Dernier, and M. D. Lukin, New J. 

Phys. 5, 113 (2003). 
[8] A. Kuklov, N. Prokofev, and B. Svistunov, Phys. Rev. Lett. 92, 



050402 (2004). 

[9] A. Isacsson, Min-Chul Cha, K. Sengupta, and S. M. Girvin, 
Phys. Rev. B 72, 184507 (2005). 

[10] A. Arguelles and L. Santos, Phys. Rev. A 75, 053613 (2007). 

[11] C. Trefzger, C. Menotti, and M. Lewenstein, Phys. Rev. Lett. 
103, 035304 (2009). 

[12] M. Iskin, Phys. Rev. A 82, 033630 (2010). 

[13] J. Catani, L. De Sarlo, G. Barontini, F. Minardi, and M. Ingus- 
cio, Phys. Rev. A 77, 01 1603(R) (2008). 

[14] G. Thalhammer, G. Barontini, L. De Sarlo, I. Catani, F. Mi- 
nardi, and M. Inguscio, Phys. Rev. Lett. 100, 210402 (2008). 

[15] B. Gadway, D. Pertot, R. Reimann, and D. Schneble, Phys. Rev. 
Lett. 105, 045303 (2010). 



